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Abstract 

The radiative corrections of order a for the charged and neutral current 
neutrino-deuterium disintegration for energies relevant to the SNO experiment 
are evaluated. Particular attention is paid to the issue of the bremsstrahlung 
detection threshold. It is shown that the radiative corrections to the total 
cross section for the charged current reaction are independent of that thresh- 
old, as they must be for consistency, and amount to a slowly decreasing func- 
tion of the neutrino energy , varying from ~ 4% at low energies to ~ 3% at 
the end of the spectrum. The differential cross section corrections, on the 
other hand, do depend on the bremsstrahlung detection threshold. Various 
choices of the threshold are discussed. It is shown that for a realistic choice 
of the threshold, and for the actual electron energy threshold of the SNO 
detector, the deduced v^. flux should be decreased by ~ 2%. The radiative 
corrections to the neutral current reaction are also evaluated. 



PACS number(s): 13.10. +q, 13.15.+g, 25.30.Pt 



I. INTRODUCTION 



Solar neutrinos from decay have been detected at the Sudbury Neutrino Observatory 
(SNO) [|1| via the charged current (CC) reaction 

Ve + d^p + p + e'. (1) 

In the next phase of the SNO experiment, currently underway, the rate of the neutral current 
deuteron disintegration (NC) 

Pe + d^p + n + Pe (2) 

will be also measured. 

From the measurement of the CC reaction rate the flux at Earth of the solar z/g was 
determined to be [0, 
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^^^o{ue) = 1.75 ± 0.07(stat)12;}?(syst) ± 0.05(theor) x lO^^cm-^s-^ . (3) 

The ®B solar neutrinos were also detected in the precision measurement by the Super- 
Kamiokande Collaboration (SK) using the elastic scattering (ES) on electrons. That 
reaction is sensitive not only to the charged current weak interaction but also to the neutral 
current interaction. From the SK measurement the *B solar neutrino flux was deduced to 
be 

$||(z/^) = 2.32 ± 0.03(stat)+°:°^(syst) x lO^cm-^g-i . (4) 

The difference between these two flux determinations, at the 3.3 a level, can be regarded as a 
'smoking gun' proof of neutrino oscillations, independent of the solar model flux calculation. 
By itself, z/gC NC scattering cannot account for the difference between @ and (p. The 
excess ES events must involve a neutrino species which contributes disproportionately to 
the NC rate. According to the oscillation hypothesis, some of the ®B solar oscillate 
into another active neutrino flavor z/^^. These u^^- neutrinos then cannot cause the charged 
current reaction, Eq.(|l|), but they can and do undergo NC scattering on electrons. Assuming 
that this is what is really happening, one arrives at the total ^B solar neutrino flux consistent 
with the Standard Solar Model . This agreement may be used as a supporting evidence 
for the oscillation hypothesis, which will be further tested by comparing the CC and NC 
reaction rates measured by the SNO experiment alone. 

The goal of the present work is the evaluation of the radiative (order a) corrections to 
the cross sections of the CC and NC reactions. Precise knowledge of these cross sections 
has obvious relevance for the determination of the ^B neutrino flux. Experimentally, one 
measures the number and energies of the electron events for the CC reaction, or the number 
of neutron events for the NC reaction, which after corrections for cuts and experimental 
efficiencies, is an integral over the incoming neutrino energies of the ^B solar neutrino flux 
(possibly modified by the neutrino oscillations) times the differential cross section. Hence 
any error in the cross section causes a corresponding error in the deduced flux. 

In analyzing the SNO CC data the theoretical cross section of Ref. was used. The 
assumed uncertainty of the calculated cross section is reflected in the theoretical uncertainty 
of the deduced flux, Eq.(^. However, radiative corrections were not applied to the CC cross 
section. 

The radiative corrections to the CC reaction (|T]) were evaluated by Towner |^. That 
analysis was recently questioned by Beacom and Parke , who noted that the total CC cross 
section for detected and undetected bremsstrahlung differ, according to the analysis of Ref. 

0. Such a difference is unphysical. The observation of Ref. [0 has lead to understandable 
confusion among experimentalists as to the appropriate radiative corrections to apply to 
the SNO data. While the published SNO result did not include any radiative corrections, 
the level of confidence in future CC/NC comparisons could depend significantly on a proper 
treatment of the radiative corrections. Thus, in what follows we revisit the analysis of Ref. 

01, in an effort to resolve the present controversy. 

While we have no quarrel with the basic treatment of the radiative corrections to the CC 
reaction in [Q, we confirm the observations of Ref. and identify the origins of the incon- 
sistency in Towner's results: (a) neglect of a strong momentum-dependence in the Gamow- 
Teller ^Si ^^Sq matrix element, and (b) improper ordering of limits involving EJ^^^ and 
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the infrared regulator. After correcting for these issues, we obtain identical total CC cross 
sections for detected and undetected bremsstrahlung. The results imply an £^j,-dependent 
correction to the total CC cross section which varies from ~ 4% to ~ 3% over the range of 
available neutrino energies. 

In addition to foregoing, we also recast the treatment in Ref. of hadronic effects in 
the radiative corrections into the language of effective field theory (EFT). Although the 
traditional treatment in Ref. ^ and EFT frameworks are equivalent, the latter provides 
a systematic approach for long-distance, hadronic effects presently uncalculable from first 
principles in QCD. As discussed in Ref. [^, matching the asymptotic and long-distance 
calculations (in EFT) involves use of a hadronic scale, Mhad, whose choice introduces a 
small theoretical uncertainty into the radiative corrections. We argue that the choice of 
Mhad made in Ref. is possibly inappropriate for the process at hand, and attempt to 
quantify the uncertainty associated with the choice of an appropriate value. Given the 
SNO experimental error, this theoretical uncertainty is unlikely to affect the interpretation 
of the CC results. It may, however, be relevant to future, more precise determinations of 
Gamow- Teller transitions in other contexts. 

Finally, for completeness, we revisit the analysis of the NC radiative correction computed 
in Ref. ||^. In this case, bremsstrahlung contributions are highly suppressed, the correction 
is governed by virtual gauge boson exchange, and the result is esssentially i?,y-independent. 
We obtain a correction to the NC cross section that is a factor of four larger than given in 
Ref. [^, which neglected the dominant graph. The implication of a complete analysis is the 
application of a ~ 1.5% correction to the tree-level NC cross section. 

Our discussion of these points is organized in the remainder of the paper as follows. In 
Section II we present our formalism for the CC radiative corrections. Given the thorough 
discussion of this formalism in Ref. [^], we restrict ourselves to only brief explanation of 
the basic formalism that is used to evaluate the corresponding FejTiman graphs and deduce 
the formulae for the differential cross section. In Section III we discuss the delicate issue of 
bremsstrahlung thresholds and the 'detector dependence' of the CC radiative corrections. 
In particular, we derive the corrections for two extreme cases of very high and very low 
thresholds, and for an intermediate, more realistic case We show in Section III where 
our results disagree with those of Ref. [§] , and trace the origin of these discrepancies. Detailed 
tabular evaluation of the modification of the differential CC cross section for the 'realistic' 
bremsstrahlung threshold is provided as well. In Section IV we consider the effects of the 
electron spectrum distortion. In particular, we consider the test of the oscillation null 
hypothesis, where the unperturbed u^, spectrum of the ®B decay is expected. In Section 
V we derive the corrections to the NC reaction rate and discuss the differences with their 
treatment in |p. We conclude in Section VI. Finally, in Appendix we collect needed formulae 
for the evaluation of the triple differential cross section (in Ee, and the angle between 
them) for an arbitrary bremsstrahlung threshold. 

II. GENERAL CONSIDERATIONS 

In the charged current neutrino disintegration of deuterons at rest in the laboratory 
frame, Eq.(|l|), the incoming v^. energy Ey, corrected for the mass difference A = Md — 2Mp = 
—0.931 MeV, is shared by the outgoing electron (energy E^), the energy of the relative 
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motion of the two protons /Mp, and by the energy of a bremsstrahlung photon (if such 
a photon is emitted), i.e., 

E, + /\ = E,+ p^Mp + (E^) . (5) 

This energy conservation condition must be, naturally, always obeyed. For the neutrino 
energies we are considering the motion of the center of mass of the protons can be neglected. 

Since radiative corrections are only few percent in magnitude, we follow Towner |^ and 
use for the 'tree level' differential cross section the formula based on the effective range 
theory (see p!0|JTl[] ) 



/ dacc \ 
[ dE, J 



tree ^ 



''^^V^WAM,PeEMHp')\\ (6) 



where for we should substitute Mp{E^ + A — Ee). It is important to remember that the 
radial integral, the overlap of the radial wave function of the two continuum protons and 
the bound state 

/(p2) = J u*^,^t{pr)udir)dr (7) 

also depends on the momentum p of the relative motion of the two protons. 

We plot in Fig. |^ the quantity |/(p^)p evaluated as in Ref. |T^, i.e., using the scattering 



length and effective range approximation as well as the Coulomb repulsion of the two final 
protons. The most important feature of the dependence is its width when expressed 
in the relevant units of p'^/Mp, the kinetic energy of the continuum protons. It is easy to 
understand the width of the curve as demonstrated in the figure. The dashed line represents 
the same |J(p^)p evaluated neglecting Coulomb repulsion as well as effective range. In that 
case a simple analytic expression obtains, 

2 (8) 

where Eb is the deuteron binding energy and the proton-proton scattering length is app = 
-7.82 fm. The value of the proportionality constant is irrelevant in the present context. Thus 
the width is determined essentially by ~ (^c)^/ (a^pMp) ~ 0.7 MeV (the term with p^ / E^Mp 
contributes very little to the width) in agreement with the more accurate evaluation. We 
explain the relevance of this width later, in Section III. 

The radiative corrections consist of two components: the exchange of virtual photons 
and Z-bosons and the emission of real bremsstrahlung photons. The Feynman graphs for 
the exchange of virtual 7 quanta and Z bosons are shown in Fig. ^. The bremsstrahlung 
graps are shown in Fig. ^. The photon emission by the moving electron is dominant ( graph 
(b) in Fig. |]), but the complete set of graphs must be considered for gauge invariance. 
The treatment of radiative corrections proceeds along the well tested lines developed for the 
treatment of beta decay (see Ref. for a review). 

Let us consider the virtual exchange corrections first. While the treatment of corrections 
involving only leptons is straightforward, those involving hadronic participants require con- 
siderable care. To that end, it is useful to adopt the framework of an EFT, valid below a 
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scale /i ~ 1 GeV. Long-distance physics (p ^ /i) associated with non-perturbative strong 
interactions are subsumed into hadronic matrix elements of appropriate hadronic operators. 
Short distance physics (p ^ /i) contributions are contained in coefficient functions C {jj) mul- 
tiplying the effective operators (see, e.g., the discussion in Ref. ||12|)- In the present case, 
the CC reaction of Eq. (|1|) is dominated by the pure Gamow- Teller transition '^Si ^^Sq. 
Thus, for the low-energy EFT, we require matrix elements of the effective, hadronic axial 
current. The resulting CC amplitude is 

M(35i ^1 So) = -^Kde7A(l - i,)uC{ii){^So\A^\^S,) + ■ ■ ■ , (9) 

where C(/i) is the short- distance coefficient function mentioned above; is an effective, 
isovector axial current operator built out of low-energy degrees of freedom {e.g., nucleon 
and pion fields); the denote contributions from higher-order effective operators; and 

where the /^-dependence of C(/i) compensates for that of the axial current matrix element, 
leading to a yU-independent result. In effect, the presence of C(/i) is needed for matching of 
the effective theory onto the full theory (QCD plus the electroweak Standard Model). 

Note that we have normalized the amplitude to the Fermi constant determined from the 
muon lifetime^, Gp = 1.16639(5) x 10"^ GeV^^ [|^. Thus, C(/i) contains the difference 



Arf^^^^-Ar, , (10) 

where Ar^*'^~'^^ contains the short distance virtual corrections to the axial vector semilep- 
tonic amplitude and Ar^ denotes the Standard Model electroweak radiative corrections to 
the muon decay amplitude. In the difference (0), all universal short-distance effects (Fig. 
2a-c) cancel, leaving only contributions from the non-universal parts of diagrams in Fig. 
2b-d. 

As a corollary, we emphasize that care must be exercised in choosing a value for the axial 
coupling constant qa used in computing (^SolA^I^Si). Typically, qa is determined from the 
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y^_G'A _ Ga{1 + Ar^) Ga, ^ a ^ v-^ /^^x 
G'y- G^l + Ar^) ^ + ^'^^ ' ^^^^ 

where Ar^ {^''^p) denotes the total radiative correction to the vector (axial vector) semilep- 
tonic amplitude. The CVC relation imphes Gy = GpVud, while the axial coupling constant 
is defined via Ga = gAGpVud- To the extent that ArJ = Ar^, the ratio A is just gA- As we 
note below, however, hadronic contributions to ArJ and Ar^ are in general not identical. 
While we speculate that the differences are considerably smaller than relevant here, arriving 
at a reasonable estimate requires a future, more systematic study. 

The asymptotic (short-distance) contributions to C(/i) have been computed in Ref. ^ 
using current algebra techniques and the short- distance operator product expansion. The 
result implies 



^This value is sometimes denoted in the literature. 
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C{fi) = 1 + 



a 
2^ 



M7 3 My 1 

3Qln-f + -ln-£ + -A(,) 



(12) 



where Q is the average charge of the quarks involved in the transition 

1 



Q = 7;iQu + Qa 



1 
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(13) 



Ag{fi) contains short- distance QCD corrections, and 6(/i) must be included to correct for 
any mismatch between the //-dependence appearing elsewhere in C{fi) and that appearing 
in the matrix element of Ax- Explicit expressions for the short-distance QCD contributions 
Ag{fJ.) may be found in Ref. 0. We note that the second term of Eq. ([121) (oc Q) arises from 
the sum of box diagrams involving (7, W) and {Z, W) pairs, while the third term arises from 
QED external leg and vertex corrections. When long-distance, 0{a) virtual effects arising 
from the matrix elements in Eq. (^) are included along with those appearing in C(/i), the 
/i-dependence of the third term in Eq. ([12|) cancels completely. 

Long-distance virtual photon contributions also contain an infrared singularity which is 
conventionally regulated by including a photon "mass" A. The resulting A-dependence is 
cancelled by corresponding A-dependence in the bremsstrahlung cross section, yielding a A- 
independent correction to the total CC cross section. In what follows, then, it is convenient 
to consider the 0{a) correction to the tree-level cross section: 



da. 



CC 



dace 



a 

l + -g 
vr 



(14) 



where the correction factor g depends on E^, and E^. as well as on E^ when bremsstrahlung 
photons are detected. This function receives contributions from C{fi), 



gl-' = 2 [C{^) - 1] 



(15) 



long-distance (p ^ yu) virtual contributions to the the axial current matrix element in Eq. 

(^), gl^^ and the bremsstrahlung differential cross section, gi,. 

In the analysis of Ref. , the long-distance contributions arising from virtual processes 
is obtained by treating the nucleon as a point-like, relativistic particle. The result is 

gl~' = ?1„ ( ^ I + 3Qln 



An f ^ 

2 Xnie 




(16) 



Here, L is the Spence function and the —3/8 is added to obtain agreement with the /3-decay 
correction and neutrino capture reaction z7g -|- p — > e"*" -|- n as calculated in Refs. |15|,|T6[. 
Note that when this —3/8 is added to A, the resulting expression agrees with the calculations 
of Refs. 
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We observe that the sum gv'" + Qv^ is independent of Mp. It does, however, contain 
the logarithm 

My 

3Ql„j^ (17) 

where Ma has been chosen in Ref. |^ as a hadronic scale associated with the long-distance 
part of the (VF,7) box diagram. Neglecting terms proportional to and me, the sum of 
the box and crossed-box diagrams depends on the antisymmetric T-product of currents 

e,uxpj d'xe''-^'S,\T[J^EM{x)rccm?Si) , (18) 

where J^^j and J^,^ denote the electromagnetic and weak charged currents, respectively, 
and where the /i and v indices are contracted with loop momentum and the lepton current. 
In order that the antisymmetric T-product appearing in Eq. (ITBp produce a Gamow- Teller 
transition, only the vector current part of J(j(j must be retained. In contrast, for pure 
Fermi transitions as considered in Ref. [§], only the axial vector charged current operator 
contributes. In that work, a choice for the hadronic scale was made based by considering 
TT /3-decay (a pure Fermi transition), and a vector meson dominance model for the axial 
vector charged current operator, leading to the appearance of the ai meson mass Ma as the 
long-distance hadronic scale. 

In the present case, such a choice appears inappropriate, since the relevant current op- 
erator is a vector, rather than axial vector current. To the extent that the vector meson 
dominance picture is as applicable to nucleons as to pions, a more reasonable choice for the 
hadronic scale would be nip. However, such a choice is unabashedly model-dependent and 
calls for some estimate of the theoretical uncertainty. On general grounds, it is certainly rea- 
sonable to choose a hadronic scale anywhere between the chiral scale A^^ = 47rF7r ~ 1.17 GeV 
and Aqcd ~ 200 MeV. Indeed, the latter choice could arise naturally from A-intermediate 
state contributions to the {W, 7) box diagrams. Thus, we replace the logarithm in Eq. (|1^ 

3^1^^ = 2-39lf ' (19) 

where the central value corresponds to Mhad = ^p', the upper value corresponds to Mhad = 
Aqcd] and the lower value is obtained with M^ad = A^. This range corresponds to a spread 
of 0.2% in predictions for the cross section Q While this uncertainty is too small to affect 
the determination of the '^B neutrino flux, it could affect more precise determinations of 
Gamow- Teller transitions for other purposes. 

The choice of Mhad amounts to use of a model for b{fi). A source of potentially larger 
theoretical uncertainties lies in possible additional, mo del- dependent contributions to this 
constant. While a complete study of these effects goes beyond the scope of the present 
work, we observe that the hadronic uncertainty cannot be finessed away using, e.g., chiral 



^We note that a similar estimate of the hadronic uncertainty in the box contributions to the Fermi 
ampUtude was made in Ref. fl^. 
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perturbation theory, since we have no independent measurements from which to fix the 
relevant low-energy constants. Moreover, the //-dependence introduced through the short- 
distance QCD correction Ag{^) must be cancelled by a corresponding //-dependent term in 
6(/t). To date, no calculation has produced such a cancellation. While the effect of this 
uncorrected mismatch between short- and long-stance effects is likely to be small, we are 
unable to quantify it at the present time. 

In contrast to the virtual corrections, the bremsstrahlung correction gi, is relatively free 
from hadronic uncertainties. In order to evaluate the bremsstrahlung part, one has to add, 
in principle, the contribution of all graphs with photon lines attached to all external charged 
particles. Only the sum of these graphs is gauge- invariant. However, for the low energies, 
the electron bremsstrahlung dominates over the proton, deuteron, and W bremsstrahlung. 

Writing again the correction to the cross section in the form l+a/ngbiEe, E^) one obtains 
the differential bremsstrahlung correction in the form 



dgb{Ee, El,, k) 
dk 



El, -|- A — Ef, — E-^ 



E, + A-E, 



k' 



2E^ 



X 



+1 



dx 



E^ 



El{E^-(3kx] 



^^,E^ + E,l 



ex^/E. 



E,, 



{E^ - (3kxy 



(20) 



{kP' + A^)^/^, i.e. A is as before the 'photon 



where k is the photon momentum, and E^ 
mass'. Also, x = cos(6'e,7). 

The dependence on the 'photon mass' A is eliminated only when one adds to the A 
dependent part of the virtual correction A an integral over the bremsstrahlung spectrum up 
to some ii^™" ^ A. We will discuss the various possible choices of E!^^^ in the next section, 
but here as an example we evaluate one of the integrals that appears in that context 



fET" k'^dk 


Jo 


VA2 + A;2 


Ej 
mg 


[h VA 


9 


'2'Cjmin 

In ^ 
A 



dx 



1 {E^-pkxf 



^min 



Pdk 



dk 



VA2 + P(A2 + ml/Elk^) 
X'^dk 



V\^ + k\\^ + ml/E^k^) 



(21) 



where f3 = Pe/Eg and the last integral, which is independent of i?™", was evaluated after 
the substition z = k/X in the limit i?™"/A oo. One must not use the limit EJ^^^ — >■ 
before all terms containing A are eliminated. 

To evaluate the full radiative correction, we assume that in an experiment one measures 
the number of events with energy Eobs ± dEots- Here Eobs = E^. when the bremsstrahlung 
photon (if such a photon is emitted) has energy of less than E'^^^. We will also assume 
that when E.^ > 



then Eobs = Eg + E^. (In the next section we will also consider 
a modification to the latter rule, making it closer to the actual conditions of the SNO 
experiment 0.) 

Thus the radiative correction to the cross section can be expressed as 



K dEobs , 



rad 



a 

TT 



'g, + gi°-{E,<E:;'- 



(22) 
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We describe in the next section how to evaluate these three functions in general and for 
three particular choices of 

III. RADIATIVE CORRECTIONS TO THE CC CROSS SECTION 

Treatment of radiative corrections involving virtual photon exchange as well as 
bremsstrahlung photon emission is a delicate issue due to the appearence of infrared di- 
vergencies. In our analysis we follow the conventional approach of introducing infrared 
regulator in the form of photon mass A, and split the bremsstrahlung contributions into two 
pieces below and above the threshold value EJ^^^ as explained above. When contribution 
from virtual photon exchange are added to the piece with < i?™", the dependence on 
infrared regulator A is eliminated. However, it is effectively replaced by a dependence on 

The threshold EJ^^^ is a detect or- dependent quantity, and may vary depending on the 
experimental conditions. In addition, the experimental conditions also dictate how to com- 
bine the piece with E^ < EJ^^^ ( + gl°^) and the E^ > EJ^^"' part. Thus it is impossible 
to give a completely general recipe here. 

With this caveat in mind, in our analysis we adopt a following framework. Each detected 
CC event is characterized by the recorded energy Eobs which in general is a function of 
electron energy E^ and, if present, photon energy Ey. Eobs = Eobs{Ee, E^). In the following 
we concentrate in particular on the role played in this context by the threshold E^^"^. 

We consider the following situations: 

A) The electrons are always recorded above the electron detection threshold E"^^^, and 
the bremsstrahlung photons are never detected, i.e. ii^™" —>■ cxd. 

B) The electrons are always recorded above the electron detection threshold and 
the bremsstrahlung photons are also always detected, i.e. iJ™" — >• 0. 

C) A more realistic case, resembling the actual situation in the SNO detector ^ when 
only part of the E^ energy is recorded, namely Eobs = {Ee — nif,)9{Ef. — i?™") + irtf. + {E^ — 
E-^i-^)e{E^ - E™"). Here e{x) is the step function. 

We simplify the cases A) and B) even further by considering an idealized detector with 
j^mm _ g^jj electrons and, thus, all neutrino interaction events are detected. When 

integrated over Eobs one arrives at the total number of events caused by a neutrino of energy 
E^. That quantity must be, naturally, independent on the bremsstrahlung threshold E!^^^. 
This is the consistency requirement imposed by Beacom and Parke [^. We verify that our 
results fulfill this condition. 

In Fig. ^ we plot as an example the normalized radiative correction to the differential 
cross-section ^"^^ / a^lf^ ( -E*^ ) for E^, = 10 MeV and two extreme cases — > oo 
(bremsstrahlung never detected, full line) and ii^™" — ^ (bremsstrahlung always detected, 
dashed line). The two corresponding curves are quite different, reflecting the different de- 
pendence of Eobs on Ee and E^. However, the areas under the curves are equal as they must 
for consistency. 

It is interesting to note that evaluation by Towner considers the same limiting cases. 
However, the results of Ref. give different corrections to the total cross section, 5(7*°*, i.e. 
they fail the consistency check. In fact, our results and Ref. differ in both extremes. We 
trace now the origin of these discrepancies. 
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A. The case of i?™" —>■ oo, no bremsstrahlung detected 



Let us first consider the limit i^;™" ^ oo. In that case we have to integrate the 
bremsstrahlung spectrum over the photon momentum from oo. At the same time, 
the energy conservation condition, Eq. (H), must be obeyed. Since now Eobs = Eg, then for 
a fixed E^, + A — Ee the quantity p'^/Mp must be varied together with E^. As noted above, 
and illustrated in Fig. |l], the quantity |/(p^/Mp)p is a rapidly varying function which falls 
off quickly for p'^/Mp > 0.7 MeV. Therefore to correctly account for this dependence we 



write 



dacc\ a (dacc\ f^'r" \I{E, + A - E, - E';)\' dg,{E,, E,, k) 



dE,J. TrVdE, /_7o \I{E, + A-E,)\^ dk 



-dk, (23) 



where gb{Ee, E^, k) is given in Eq. (|20D. If \I{Ej, + A — E^ — Ey)\'^ could in fact be treated 
as a constant, the ratio of the two P would be unity, and Eq. (|^) would be identical to Eq. 
(13) in 0. To make the connection with Ref. ^ even more concrete we write (note that for 
E < K™" R , = E )■ 



da, 



cc 



dE, 



obs / h \ ^-^obs . 




b 



f^- (im + A-Kb^-E^ \ dgbiE^bs,E^,k) 
Jo [ \m + A-E,bsW ~ ) dk ■ ^^^^ 

(In Eqs. ( [23|J2^ ) the upper limit of the integral obviously should not extend beyond the 
corresponding bremssstrahlung endpoint.) 

The first term on the r.h.s of Eq. (p^ is the contribution present in Ref. 0. It contains 
infrared divergence that disappears after contributions from virtual photons are added. The 
second term is infrared finite. Due to the shape of \I{Ei, + A — E^. — E!^)\'^ this term 
enhances the contribution of the low energy tail in {dacc/dEobs)b- The overall result of the 
low Eobs tail enhancement is that the total cross-section is increased by about 3% compared 
to corresponding result in for the considered case of E^, = 10 MeV. 

B. The case — > 0, bremsstrahlung always detected 

If one wants to study the opposite extreme EJ^^^ ^ it is crucial in Eq. (p^) to first add 
all three terms, eliminate infrared cutoff dependence, and only then take the limit E"!^^^ —>■ 0. 
The order of limits A — ^ and is important because upper limit of integrals like Eq. 



(pi|) is -E™"/ A. Since A ultimately is an infinitesimal unphysical parameter it is mandatory 
to maintain EJ^^^ 3> A during the entire course of the calculation 0. This leads to a non- 
intuitive result that the second term in Eq. (^2]) has a non-zero contribution even in the 
limit i^™" — i> 0. In particular, one must write the following expression corresponding to the 
second term on the r.h.s of Eq. (p^), i.e. for E^ < i?™": 



3lf A were truly the photon mass, the requirement that > X would be obvious. 
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a 



2 In 



tree 



TT 



1 , 



(25) 



C(/5) = 21n(2) 

■1-/52 



2/? \l-(3 



4/3 I 1-/3^ 



2 + In 



4 



+ i[L(/3)-LM)l + i5(L(l^^ 



'1 + /5' 



where for (3 <1 the Spence function is 

f/3 ln(|l 



00 ok 



dx 



X 



E 

fc=i 



A;2 



(26) 



(27) 



The A-dependent terms in Eq. (|25f) will be cancelled by A-dependent pieces from virtual 
photon contributions, and the logarithmic divergence in E^^^ will disappear after the piece 
with > E!^^" is added to the cross-section (third term in Eq. (|22|) ). Only after this is 
done is one allowed to take EJ^^^ — > 0. The most striking feature of Eq. ( p6D is that it 
is independent of E"^^^. Consequently, it survives in the limit E"^*"- — ;> 0. It appears that 
this procedure was not followed in Ref. and, therefore, equation (44) and Table II in |^ 
must be modified accordingly. We plot in Fig. ^ the cross section correction {a/TT)g{Eobs), 
which is for — > independent of the neutrino energy E,y. Note that it differs in slope 
compared with its analog in Table II of Ref. . 

The two aforementioned modifications to the treatment in Ref. ^ allowed us to bring 
the two cases E"!^^^ —>■ 00 and ii^™" — > in agreement in terms of correction to the total 
cross-section and resolve the discrepancy mentioned in Ref. 0. In either of these extreme 
cases, by integrating over Eobs we obtain the QED correction to the total cross section as a 
function of the neutrino energy E^,, 5a^°^{Ey), displayed in Fig. ||. 



C. Realistic bremsstrahlung threshold 

The treatment of the more realistic case is now straightforward. First and second terms 
(virtual and < i^™") on the r.h.s of Eq. ( ^2]) are evaluated by setting i^;™" = 1 MeV 
and Ef. = Eobs- In the third term one has to set E^ + E^ = Eobs + i?™" = const and 
integrate over Eg. 

In particular, suppose we write the double differential cross-section for the + z/g — >■ 
p + p + e + 'yas d^acc / {dEedE^) = f{Ee, E^). Then the total cross-section with E^ > i^;™" 
MeV is: 



{alc\ot= / dEj f{E,,E,)dE, 

J me J EJ^"-" 

dEobs I fiEe, Eobs + ^7" - E,)dE, . (28) 



We simply performed the change of integration variables from [Eg, E^) to [Eobs, E^) in 
the spirit of Ref. 0. Now we can write in the notation of Eq. (^21): 
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a 



high^^^ > ") = r"f{E,, E,,s + - E,)dE, . (29) 



The result, as expected, is a function of Eobs only. In order to generalize to the case 
rrie one has to be careful because change of variables from {Ee,E^) to {Eobs, Eg) 
becomes less trivial. It is possible to show, however, that the following relationship holds: 



-Qb 1-^7 > E > me) - —gi, [Ey > h - m^) 

71 ' 71 ' 



J^rmn 

+ / ^ [f{E,, Eobs + ^7" - me) - f{E,, Eobs + - E,)]dE, , (30) 



where f{x,y) is the function defined before in Eq. (p8|). Eq. (^) allows one to obtain the 
correct spectrum, Eq. (P^, for any electron threshold in terms of the ideal case where all 
electrons are detected. We note that it is only the third term in Eq. ( p2[ ) that (implicitly) 
depends on The impact of the refinement in Eq. (0) is rather small for low values of 

E™". We evaluated it for E^^"- = 1.5 MeV (1 MeV kinetic energy). The effect of the second 
line in Eq. (^Oj) is a 0.03% modification of the differential cross section. Consequently, we 
neglect this refinement in our analysis. 

The spectrum for the case C) is shown in Fig. ^ in the dash-dotted line. As expected, 
the upper 1 MeV of that spectrum coincides with the i^™" 00 case. Note that the areas 
under all three cases in Fig. ^ are the same, as they must for consistency. 

In Table | we provide detailed tabular information on the correction to the differential 
cross section for the full range of neutrino energies E^, and Eobs for the case C). 



IV. FOLDING WITH THE '^B SPECTRUM. 

In an actual solar neutrino experiment, like SNO, the recorded quantity is the number 
of events with energy Eobs (or the total number of events integrated over Eobs)- This is 
an integral over the product of the incoming neutrino spectrum and the differential cross 



section. I.e. 



dEobs Jo dEobs 

where f{Ey) is properly normalized incoming neutrino spectrum, possibly modified by neu- 
trino oscillation. When testing the 'null hypothesis', i.e. asking whether neutrinos oscillate, 
one takes for the incoming neutrino spectrum simply the shape of the v^, spectrum from 
decay [|l7| (normalized to unity over the whole range of E^). 

In Fig. we show the folded correction to the differential cross section, Eq. (|3l|) (full 
line). The case C), i.e. the realistic bremsstrahlung detection threshold has been used to 
produce the full curve. For comparison we also show similarly folded tree level cross section, 
scaled by a factor 1/40 so that it fits in the same figure (dashed line) . One can see that the 
two curves are similar in shape which is basically dictated by the incoming spectrum, 
but the QED correction is shifted toward smaller Eobs-, roughly by the value E"^^"^ = 1 MeV. 

When integrated from the threshold used in the SNO analysis, -EX" ~ '"^e = 6.75 MeV, 
the full line represents ~ 2% increase of the total total cross section, therefore ~ 2% decrease 
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of the deduced flux, Eq. (j^), when the radiative corrections are properly included. If it would 
be possible to reduce the threshold to very low values, the reduction of the flux would be ~ 
3%. 

These relative increases of the total cross section obviously differ somewhat from the 
values displayed in Fig. ^ that were obtained for monochromatic neutrinos. The difference 
is, naturally, caused by the effect of the shape of the radiative correction to the differential 
cross section in combination with the shape of the z/g spectrum. In particular, for the 
actual SNO -EX" threshold one could have expected an increase of the cross section (or 
count rate) due to radiative correction of ~ 3% based on Fig. |^ while the folding with the 
incoming spectrum reduces this value to ~ 2%. 



V. RADIATIVE CORRECTIONS TO THE NC CROSS SECTION 



The NC cross section is governed by the effective four fermion low-energy Lagrangian 



£V—had 

where 



2^2 



^T=1t/T=1 



v: 



A 



T=l 



T=l _ ^ 
1 



U7^u - (i7^(i 



1 r 



M 2 
4^=° - - 



M7^M + d-ff,d 



M7/.75M + c?7^75C? 



(32) 

(33) 
(34) 



and where only the effects of up- and down-quarks have been included. At tree level in the 
Standard Model, one has 



2(1-2 sin^ Ow) =° = -4 sin^ 9 



w 



T=0 







(35) 
(36) 



As noted in Ref. 0, the incident and scattered neutrinos do not contribute to the 
bremsstrahlung cross section at 0{G^pa), while radiation of real photons from the par- 
ticipating hadrons is negligible. Thus, the dominant radiative corrections involve virtual 



exchanges, which modify the C,va from their tree-level values: 



(37) 



where the RyA contain the 0{a) corrections. Since the NC amplitudes are squared in 
arriving at the cross section, the total correction to the NC cross section will go as twice the 
relevant RyA- 0-^ notation of Ref. R^^^ = {a/2'7i)g^'~^ .) 

As emphasized in Ref. 0, considerable simpliflcation follows when one considers only 
the dominant break-up channel: ^Si(T = 0) ^^So(T = 1). As a AT = 1, pure spin- 
flip transition, this amplitude is dominated at low-energies by the Gamow- Teller operator. 
Magnetic contributions are of recoil order and, thus, v/c suppressed. Consequently, we need 
retain only the Aj^^^ term in Eq. (^) and consider only the correction R^^^- 
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The source of corrections to R^^^ include corrections to the W and Z-boson propagators 
(Fig. ^ graph (a)); electroweak and QED vertex corrections to the Zuu and Zqq couphngs 
(Fig.^, graph (b)); external leg corrections (Fig. ^ graph (c)); and box diagrams involving 
the exchange of two W^s or two Z's (Fig. ||, graph (d)). The presence of VT-boson propagator 
corrections arises when the NC amplitude is normalized to the Fermi constant Gp determined 
from muon decay. Only the difference between the gauge boson propagator corrections enters 
the RvA tliis case. Note that Z-7 mixing does not contribute to R^^^ since the neutrino 
has no EM charge and the photon has no axial coupling to quarks at = 0. Similarly, one 
encounters no Z7 box diagrams for neutrino-hadron scattering. 

In the analysis of Ref. |Q, only the ZZ box contribution was included, yielding a cor- 
rection R^^^ ~ 0.002. Inclusion of all diagrams, however, produces a substantially larger 
correction. From the up-dated tabulation of effective u-q couplings given in Ref. [pT^ 
obtain 



we 



R 



.T=l 



PuN + ^dR - Kr + Kl - AdL - 1 ~ 0.0077 



(38) 



where we have followed the notation of Ref. |T^. In particular, the WW box graph con- 
tributes roughly 80% of the total: 



RT {WW -box) 



5a 



877 sin 9w 



0.0063 . 



(39) 



The net effect of the total correction is therefore g^'" = 6.63, i.e., ~ 1.5% increase in the 
NC cross section, as compared to the 0.4% increase quoted in Ref. H. 



VI. CONCLUSIONS 

The radiative correction of order a for the charged and neutral current neutrino- 
deuterium disintegration and energies relevant to the SNO experiment are consistently eval- 
uated. For the CC reaction the contribution of the virtual 7 and Z exchange is divided 
into high and low momentum parts, and the dependence on the corresponding scale ~ 1 
GeV separating the two regimes is discussed in detail. For the bremsstrahlung emission we 
discuss the important role of the bremsstrahlung detection threshold EJ^^^. In particular, 
we consider the two extreme cases EJ^^"' 00 and i^^*" — > as well as a more realistic 
intermediate case. We show that our treatment, unlike Ref. gives consistent, i.e. iJ™"-- 
independent, correction to the total cross section, shown in Fig. ^. This correction, slowly 
decreasing with increasing neutrino energy E^, amounts to ~ 4% at low energies and ~ 3% 
at the end of the '^B spectrum. 

The magnitude of this correction is in accord with the correction to the inverse neutron 



beta decay, + p n + evaluated in Refs. |15,16| and with the correction for the 



pp fusion reaction evaluated in Ref. [19|. Note that in these references only the 'outer 



radiative corrections', i.e. only the low momentum part of the virtual photon exchange was 
considered. The high momentum part, which is independent of the incoming or outgoing 
lepton energies, and which is universal for all semileptonic weak reactions involving d ^ u 
quark transformation, amounts to ~ 2.4% pO| and should be added to the results quoted 



in Refs. 
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We identify the origin of the inconsistency of the treatment of Ref. 0: (a) neglect of a 
strong momentum-dependence in the Gamow- Teller ^Si — ^-^Sq matrix element, which affects 
the case of i^™" oo, and (b) improper ordering of limits involving and the infrared 

regulator, which affects the case of 0. For the more realistic choice of i^^*"- we 

provide a detailed evaluation of the correction to the differential cross section. 

We also discuss the effect of folding the cross section with the (unobserved directly) 
spectrum of the ®B decay. We conclude that for the realistic choice of EJ^^"^ and for the 
electron detection threshold of the SNO collaboration, the solar Uf. flux deduced neglecting 
the radiative correction would be overestimated by ~2 %. 

Next we consider the effect of radiative corrections to the neutral current deuteron dis- 
integration, so far not analyzed by the SNO collaboration. In that case the radiative correc- 
tions, associated with the Feynman graphs in Fig. |], are dominated by the virtual Z and W 
exchange, in particular be the box graph in Fig. |^ (d). The corresponding neutrino energy 
independent correction to the NC total cross section is ~ 1.5%. 

We provide in the Appendix a set of formulae relevant for the case of an arbitrary 
bremsstrahlung threshold EJ^^^. These formulae allow one to evaluate the CC reaction 
differential cross section in terms of the electron energy Ee, the photon energy E^, and the 
angle between the momenta of the electron and photon. 
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Appendix 

Here we provide a recipe for obtaining radiative correction to differential cross-section 
for the reaction d + i>e^e + p + p. The prescription is infrared finite and allows arbitrary 
values of cutoffs for detection of electrons and photons. 

Unlike the discussion in the main text here we do not choose any particular model 
for what an experiment can detect. The only assumption that we make is that the 
bremsstrahlung photons cannot be seen below certain energy EJ^^"^. Therefore, contribu- 
tions from all photons with energies below this cutoff are added and the only thing that is 
left for detection is the electron energy. 

We make no assumptions as to how photons with > are recorded. For their 

contribution we provide triple differential cross-section that depends on electron energy, 
photon energy and angle between the direction of the electron and emitted photon. This 
expression can be incorporated in the detector-specific simulation software for appropriate 
analysis. 

We combine the contributions from photons with E^ < E^^^ with virtual photon and Z 
exchanges to get infrared finite answer (first two terms in Eq. (|2^)). The result is: 



/ da{Ee,E,) 
V dE, 



[ dE, 



tree 



-{2 In 



E„ 



nip 



+h{Emax, E,) + h{E,ma., E,) + C{(3) + A {(3) 



Er, 



+gv + -In — + 3Q\n 



E^a. f\I{E^ + A- E, 



obs 



E^ 



V \I{E, + A - Eobs) 



A dgbjEpbs, Eu, E^) 



dE^ 



Min 



e: 



E^ + A-E, 



(40) 
(41) 



where C(/3) is defined in Eq. ([26|), gv^ is taken from Eq. (p!5|), and A'{l3), Ii{E., 
and l2{Emax, Ey) are defined as follows (see 0): 



E ) 

max ) } ■) 



HE, 



maxi Ey^ 



h{E^, 



E ) 



P 



1 



-In 



PEl 



In 



' 1 + I3 \ {E, + A- E, 
I- (3) 15 



5-3 1 



Er, 



E„ + A- Ep 



1 





In 



'1 + /3' 
1 



- 1 



Q 



E^ + A-E, 



E„ 



E„ + A- Ep 



/ O 



3/2 



- 1 



dQ 

Q 



(42) 



For E^ > EJ^^^ we write the triple differential cross-section 
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dEg dE^ ^•^/ 



TT TT 



x{Mp{E, + A-E,- E^)y/^\I{E, + A-Ee-E^ 
1 



xE^ 



E^,{l-Px) 



2 Eg + Ej 1 



(43) 



where a; is cosine of the angle between photon and electron momenta. We have integrated 
over the corresponding azimuthal angle. 



17 



TABLES 



TABLE L Values of the correction to the differential cross-section, 6 dE\°'' ) /^^"^i^'^)^ 
normalized to the total tree-level cross-section, in %/MeV. The neutrino energy Ei, in MeV labels 
the columns, while the total energy observed in the detector, Eobs, in the form E,^ + A — Eobs also 
in MeV, labels the lines. The dash-dotted curve in Fig. ^ corresponds to the column with £"2^=10 
MeV. 
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FIGURES 




p7Mp (MeV) 



FIG. 1. Radial integral |/(p^/Mp)p. Exact result (solid line) and scaled approximation with 
Coulomb repulsion and effective range set to zero (dashed line) are shown. 
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(c) (d) 

FIG. 2. Order a radiative eorrections to the charged-current breakup of deuteron: 
d + fe p + p + e involving virtual 7 quanta and Z boson exchange. The large shaded oval 
represents the vertex with all its hadronic complications. 
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(b) 

FIG. 3. Order a corrections due to bremsstrahlung emission. See caption to Fig. |2|. 
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FIG. 4. Corrections to the differential cross-sections as a function of the observed energy, 
normahzed to the total tree-level cross-section. Solid line corresponds to E'^'^'^ — >■ oo (-Bobs = E^), 
dashed line to EJ^^^ — (£?o6s = E^, + Ej), and the dot-dashed line is obtained by setting 
Eobs = Ee + {E^ - l)9{Ej - 1). All lines are evaluated for E„ = 10 MeV. 
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FIG. 5. Correction to the cross section for — ^ 0. Note that in this case the corrections 
does not depend on Ei,. 
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FIG. 6. Radiative corrections to the CC total cross section as a function of neutrino energy. 



5.50 




FIG. 7. Radiative corrections to the total CC cross section folded with the incoming °B Vf. 
spectrum as a function of the detected energy Eo\,g (full line) Also shown is similarly folded tree 
level CC reaction cross section, scaled by a factor of 1/40 (dashed line). 
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(d) 

FIG. 8. Feynman graphs relevant for the radiative correction to the NC cross section (see text 
for explanation). 
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